Chiral entanglement in triangular lattice models 
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We consider the low energy spectrum of spin-j two-dimensional triangular lattice models subject 
to a ferromagnetic Heisenberg interaction and a three spin chiral interaction of variable strength. 
Initially, we consider quasi-one dimensional ladder systems of various geometries. Analytical results 
are derived that yield the behavior of the ground states, their energies and the transition points. 
The entanglement properties of the ground state of these models are examined and we find that the 
entanglement depends on the lattice geometry due to frustration effects. To this end, the chirality of 
a given quantum state is used as a witness of tripartite entanglement. Finally, the two dimensional 
model is investigated numerically by means of exact diagonalization and indications are presented 
that the low energy sector is a chiral spin liquid. 
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I. INTRODUCTION 



Exotic quantum orders of correlated many-body sys- 
tems have been studied extensively over the lastyears [l| . 
Prime examples include the Laughlin states [2|, which 
were studied in relation to the fractional quantum Hall 
effect, and the chiral spin states Q that were studied in 
relation to high-T c superconductivity. Chiral spin states 
break spontaneously the parity (P) and time-reversal (T) 
invariance. It is known that in two spatial dimensions it 
is possible for such "incompressible" quantum systems 
with a mass gap to exhibit topological order One 
of the most intriguing properties in this case is that the 
degree of the ground state degeneracy depends on the 
topology of the surface in which the system resides. As 
a result, topologically ordered quantum states cannot be 
distinguished by any local measurements but are only 
globally distinct. The outstanding problem of character- 
izing quantum states that possess topological order [5| 
has recently motivated the study of the many-particle 
entanglement properties of such states @. 

In this paper we consider a model that breaks the PT- 
symmetry explicitly via a scalar chiral interaction term 



X = a 



i ■ <7j X <Jk 



(1) 



while it preserves the SU(2) rotational symmetry. We 
study the quantum phase transitions of this model and 
the properties of its different phases. In particular, we 
are interested in the local indistinguishability of the de- 
generate ground states, and the absence of any structure 
in the classical correlations jx, Q ■ In addition, we intro- 
duce an entanglement witness that quantifies the quan- 
tum correlations of the system due to the chiral currents. 
This is a suitable multipartite entanglement witness that 
detects two-body and different classes of three-body en- 
tanglement. 



From the experimental point of view, measurements 
of the scalar chirality in pyrochlore ferromagncts 8] 
and Kagome lattice structures @ have been performed. 
Moreover, recent theoretical and experimental advances 
in cold atomic and molecular physics indicate the pos- 
sibility of generating the required chiral interactions. 
As we shall see, the chiral model studied here could 
be implemented with cold atoms superposed by opti- 
cal lattices in the presence of effective magnetic fields 

nam Em, p. 

The model we study is a two dimensional triangular 
lattice of spin- ^ particles subject to chiral and ferromag- 
netic Heisenberg interactions. The Hamiltonian of the 
system is given by 



H(\) 



(id) 



Sj + A 



E 



(2) 



where and (i,j,k) denote any two and three near- 
est neighboring sites, respectively, and $i = (Xi,Yi, Z{) 
where Xi,Yi,Zi are the Pauli matrices for spin i. The 
parameter A is a real number that determines the rela- 
tive strength of the chiral and ferromagnetic interactions. 
The sign of A is irrelevant for the eigenvalues and it can 
be changed with a time reversal transformation. 

The ferromagnetic and chiral terms favor different or- 
ders, although they are both SU(2) symmetric. When 
|A| <C 1, the Hamiltonian @ has a ferromagnetic ground 
state which is a product state with all spins pointing 
along the same direction. On the contrary, when |A| ^> 1, 
the chiral part favors a ground state which is entangled, 
typically belongs to the singlet sector and has zero mag- 
netization and nonzero chirality. 

The focus of our work is on the properties of the chiral 
phase and on the phase transition that takes place as the 
value of |A| is changed. Our study begins in Sect. HT1 with 
an interpretation of the chirality for three spin systems, 
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TABLE I: Half of the eigenstates of the chirality operator of 
Eq. {1} for three spins expressed in terms of the phase to — 
exp(27ri/3) and its complex conjugate uo* . The eigenstates 
with negative values of S z are obtained by nipping all spins. 

focusing on the relation between chirality and entangle- 
ment. In Sect. Unl we introduce simple quasi-lD geome- 
tries, such as regular polygons and spin ladders, and ana- 
lyze the ground state both numerically and analytically. 
We present indications of a topological quantum phase 
transition and study how different levels of frustration 
on the chirality give rise to ground states with different 
kinds of entanglement. In Sect. IIVI we move to a more 
sophisticated geometry, namely a 2D lattice with peri- 
odic boundary conditions which forms a torus. Here we 
perform a numerical study of the phase transition and 
characterize what we conjecture to be a spin-liquid or- 
der in the chiral phase. Finally, as commented above, 
Sect. [V] discusses the possibility of realizing this model 
using cold atoms in optical lattices, and we summarize 
our conclusions in Sect. [Vlj 



II. CHIRALITY AND ENTANGLEMENT 

Consider a system of three spins in a triangular con- 
figuration. We pose the following question: is there any 
relation between entanglement properties of a pure state 
of these spins, \^), and their scalar chirality, 

X = <*l*l*>, (3) 

defined as the expected value of the operator in Eq. (|TJ) ? 
In subsections B and C below we answer this question 
and show that indeed the chirality can be used to detect 
the existence of two- and three-partite entanglement in 
arbitrary states. This is true for any number of spins, 
which means that by computing the chirality on differ- 
ent sets of spins one can learn about the distribution of 
entanglement in the ground state of a Hamiltonian, such 
as the one of Eq. ([2]) (cf. Sect. IIII[) . Clearly, the scalar 
chirality can also be defined for general mixed states, p, 
as tr(pX). This definition is employed in Sections IIIII 
and IIVI to characterize the multiqubit ground states of 
quasi- ID or 2D systems. 

A. Eigenstates and eigenvalues 

The chirality operator X is Hermitian, imaginary - 
i. e. changes sign under complex conjugation — and it is 



invariant under global rotations. The Hermitian nature 
implies that \ is real, while the rotational symmetry im- 
plies that we can split its eigenstates into the eigenstates 
ofS 2 andS z . 

From the relation A" 2 = — (<7i+<7 2 +<73) 2 + 15 we deduce 
that the chirality operator can only take values between 
+2\/3 and —2^/3. This relation is further exploited in 
Ref. [3j to compute all eigenstates and eigenvalues, which 
we summarize in Table HI The eigenspace of X decom- 
poses into a spin- 1 multiplet with eigenvalue x = 0, and 
two different spin-^ with eigenvalues \ = ±2-^/3, so that 
in total we have 4 + 2 + 2 = 8 eigenstates, as expected. 

Remarkably, the states with nonzero chirality are ob- 
tained when a single spin is flipped with a relative phase 
lo or to* between sites. Such a wavefunction can be inter- 
preted as a current moving left- or rightwards thus hon- 
oring the name "chirality" of the operator. This point 
is even more evident in the bosonic model developed in 
Sect. IIII Bl but it is also a signature of the relation be- 
tween chirality and W-like entangled states (l5j . 

B. Three-partite states and chirality 

The space of states for three spins (or qubits) is a par- 
ticular case in which all states can be classified into few 
classes depending on their entanglement properties [l5l |. 
In this section we construct explicitly the most relevant 
instances of each class and analyze the corresponding val- 
ues of the chirality. The goal is to find a relation between 
chirality and pure three-partite entanglement. 

To begin with, we consider product states of the form 
I*) = |t) (8 (o|t> + b\l)) 8 (o'|t> + b'\l)), where without 
loss of generality the first spin is oriented upwards. We 
find that x = — 4Im(a*6a'6'*) becomes maximum when 
\a\ = \b\ = \a'\ = \b'\ = l/v2 and the phases are suitably 
chosen. Thus, for a product state we will always have 
Ixl < 1- If the chirality takes a value larger than I, it 
signals the presence of quantum correlations between two 
or more spins. 

We now consider bipartite entangled states of the form 
| ^ 12 ) <8> | ^3 ) • Due to rotational symmetry we may assume 
|* 3 ) = |t> and |M> 12 ) = a |tt)+&|t4)+c|lt)+d|U). Con- 
sequently, we obtain x — 2i(bc* — b*c), with a maximum 
absolute value given by |x| = 4|6c| < 2. It is noted that, 
in this case, |y| is equal to twice the concurrence of the 
I M/12) state [16}. To see that this is not accidental con- 
sider the state 1*12) = «|tt) +/3|4D- The local unitaries 
that optimize the expectation value of the relevant part 
of the chiral operator, X1Y2 — Y1X2, implement the ro- 
tations X2 — > — Y2 and Y2 — > X2. Thus, the maximum 
expectation value of the chiral operator corresponds to 
the expectation value of 2Y1Y2, which is twice the con- 
currence. Continuing with our reasoning, we have found 
that the maximum chirality bipartite entangled states 
can give is x = 2. If the scalar chirality assumes a larger 
value, then genuine three-party entanglement must be 
present. 



We finally arrive to the three-partite entangled states, 
which can either be GHZ-like or W-like states [l5j . In 
the case of the GHZ-like state |*) = a|tTT> + b\ili), we 
obtain \\\ — |a6|3-\/3 with maximal value x — 3v3/2 > 2. 
In the case of W-like states the maximum value of x is 
attained for the eigenstate of the chirality operator with 
the maximal eigenvalue, \\,\) , which gives \ — 2\/3 
(see Table H]). Thus, we can employ the chirality of three 
spins as an observable to distinguish between the different 
types of three partite entanglement. 

It is noted that for states which are symmetric under 
the exchange of any two spins, 1, 2 or 3, the above criteria 
simplify. Namely, it can be proven that the chirality x 
is different from zero only if the state possesses genuine 
tripartite entanglement. 



C. Chirality operator as an entanglement witness 

We showed in the previous subsection that for all sep- 
arable states I*) we have |(*|-^|*)| = |x| < 1. States 
whose chirality is greater than 1 necessarily have some 
amount of entanglement. Thus the chiral operator X can 
serve as a witness of entanglement [I?} • An operator O is 
a suitable entanglement witness if for all separable states 
|\&) its expectation value (O) = (^\0\^) is bounded and 
only entangled states can yield expectation values that 
exceed this bound. 

In our case, the appropriate definition of the chiral 
entanglement witness is 

E x (i&) = -l + max|xi >j , fc | (4) 

where the maximization is over all unitary operators 
U G SU{2) acting locally on each of the spins k. If 
Ex{^) > then |^) is guaranteed to be entangled, other- 
wise we cannot infer the presence of entanglement on this 
basis alone. It should be stressed that Ex(^) > 1 implies 
the existence of tripartite entanglement, and that, in con- 
trast to previous other entanglement witnesses [181 ] , it can 
distinguish between some of the GHZ- and W-entangled 
states, as shown in the previous subsection. Finally, note 
that the definition of E x can be straightforwardly gener- 
alized to mixed states. 



III. QUASI-ONE DIMENSIONAL GEOMETRIES 

In this section we consider quasi-one dimensional (ID) 
geometries of ladders and polygons as a first step towards 
studying the two dimensional (2D) case. We will estimate 
the ground state properties of Hamiltonian ^ using a 
mean field variational estimate and verify that it agrees 
with the outcome of exact numerical diagonalizations. 
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FIG. 1: Quasi-ID geometries. The spins reside on the vertices 
and spin-spin interactions are represented by the edges. At 
each triangle there is also a chiral three-spin interaction which 
has the same sense everywhere, as indicated by the arrows. 

A. Different types of quasi-lD spin systems 

We will particularize the model of Eq. |J2) to four differ- 
ent geometries which include three types of spin ladders 
and a ring-type setup (see Fig. QJ. It is useful to repa- 
rameterize the Hamiltonian as follows 

H = — Jj,j(7i ■ <7j + A ^i,j,k^i ' <?j x <?fc- (5) 

i<j i<j<k 

Here Jij and -X"j,j,fc are nonzero only for the appropriate 
bonds of each particular lattice and determine both the 
geometry of the lattice and the sign of the Heisenberg 
and chiral terms ( Jij , X^j^ G {0, +1, — 1}). 

For |A| <C 1, the ferromagnetic interaction dominates 
and spins align along any direction. For |A| ^> 1, the 
chiral interaction dominates, the magnetic ordering dis- 
appears and we observe the establishment of chiral order 
in different plaquettes. Note that by fixing the sign of 
Xijk we impose the condition that all plaquettes have 
the same chiral orientation [See Fig. [T] . This leads to an 
additional frustration in the chiral regime and may pre- 
vent the system from achieving the largest value of the 
chirality, which is 2^3 per plaquette, even for very large 
values of A. 

Let us now look at the different quasi-lD models. For 
the ladder of type A the nonzero couplings are 

J M+ i = 1, (6) 

JiA+2 = 1. 
-X*i,i+l,i+2 = (—1)', 

where the factor (— l) 1 fixes the same orientation of the 
chiral term on all triangles. The ladder of type B is very 
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similar, but the upper row of Heisenberg and half of the 
chiral interactions are missing, so that 



Ji,i+1 = 1) 
J2i,2i+3 = 1) 
-^2i,2i+l,2i+2 = 1- 



(7) 



By contrast, the ladder of type C is formed by a set of 
weakly connected triangles, in which case we obtain 

J M+1 = 1, i = l,2,4,5... (8) 
Ji,i+2 = 1) i = 1,4, 7, ... 

JiS+3, = 1) 4 = 2,5,8,... 
^Q,i+l,i+2 = 1) 4 = 1,4,7,... 

Finally, in the ring-type geometry shown in Fig. [TJD) 
there is a central spin-i particle (the spin labelled '0') 
which is connected to N— 1 equidistant spins. The Hamil- 
tonian of this system is best written explicitly as 



-Bring — — <?0 ' Oi 



N—l 

»=i 



iV-l 

+ A ^ (To • Si X (T. i+ l. 



(9) 



B. Analytic results 

We have studied analytically these models, in order 
to get information about the ground state energies and 
the location of possible quantum phase transitions. We 
present in detail the method only for the type-A ladder, 
but the procedure is the same for the rest of the geome- 
tries. Our solution works in three steps (i) a mapping 
to a hard-core bosonic problem, (ii) a second mapping to 
a fermionic model and (iii) a mean-field solution of the 
fermionic model. 

First of all, the spins are mapped onto bosons using 
the Dyson-Maleev transformation 

X n = b n + bl Y n = -i(b n -bl), Z n = 1 - 2&t&n (10) 

where the b, b^ are bosonic operators that satisfy the 
commutation relations = S nm , and [b n ,b m ] — 

[6jj,fcj„] =0. To be consistent with the spin-i problem 
at hand, we allow up to one boson per lattice site, i.e., 
tfibn < 1. In this representation the Heisenberg interac- 
tion between the m-th and ro-th spins reads 



(1 - 2blb m )(l - 2&t bn ) + 2 (blb n + b\b m ) 



(11) 

while the chiral interaction between spins Z, m and n 
becomes 

X lmn = 2i[{blh - b\b m ){l - 2blb n ) 
+(bib m -bl n b n )(l-2bjb l ) 
+ (blb n -bib l )(l~2blb m )}. (12) 



The next step is to turn to fermionic variables, c„. We 
achieve this by using the Jordan- Wigner transformation 



Cn = (-l) C "kn, Cn = ^2 



C' k C k 



(13) 



k<n 



where £ n gives the total number of fermions on the 
sites to the left of site n. These new operators sat- 
isfy fermionic anticommutation relations {c n , c^} = 5 nm , 
and {c n ,c m } = {cj, , c^} = 0. In addition to this we also 
need the following relations 



bib^^cU-ir^^c,^. 



(14) 



Replacing all formulas above into the particular Hamil- 
tonian ([3]) one obtains a complicated fermionic model 
which in general has no simple analytical solution. We 
may nevertheless estimate variationally the properties of 
the ground state. 

Our approximation method proceeds as follows. We 
notice that for |A| 1 one of the possible ferromagnetic 
ground states is polarized along the z-direction. This 
means that the state has no effective fermions and a = 
(b^bn) = 0. It thus makes sense to treat a as our order 
parameter; to accurately find the phase transition the 
expression as a function of a must become exact in the 
limit of one boson, a = 0(1/N), where the transition to 
the chiral regime occurs. 

Our mean field method dictates two approximations. 
In the Heisenberg term we replace 

(l-2c, t 1 c n )(l-2c I t „c m ) -> l+2(a-l)(4c„+4c ro ) (15) 

for the various pairs of interacting sites. In the chiral 
term we simply substitute {2c\ +2 c i+2 — 2cj_ 1 Cj_i) with 
0, which leads to a cancellation of the diagonal hopping 
terms, leaving only the quadratic next-nearest-neighbor 
hopping. With this, the Hamiltonian of type-A ladder 
becomes 

H A = -2N -2^2 [4+i c l + c i c '+i ~ 8 (« - l ) c \ c i 
i 

-2 {[(1 - 2a) + A4(-l)'h f +2Ci + H.c.} . (16) 
i 

It is now convenient to regroup the fermionic operators 
to form a spinor, = (c2/+i,c/), and then to perform 
a discrete Fourier transform 



(17) 



where p = 2-KnfL for n = 0, 1, . . . (L— 1) and L = N/2 is 
the length of the ladder. Consequently the Hamiltonian 
becomes 



H A = -2N + J2 K M p^ 



(18) 
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with the coupling matrix 

M p = [8(1 - a) - 4(1 - 2a) cos(p)]I 

— 2[1 + cos(p)]<Ta; — 2sin(p)(T y — 4Asin(p)er 2 . (19) 

The eigenvalues of the matrix yield the eigenenergies of 
the Hamiltonian, which are 

E± (?) = 8(1 - a) - 4(1 - 2a) cos(p) 

±2^/(4A 2 + l)sin 2 (p) + [1 + cos(p)] 2 . (20) 

If we label as the 2L eigenenergies, sorted such that 
Ek < Ek+i, the lowest energy state of our model has aN 
fermions occupying the states k — 1 up to k = aN. The 
energy is then given by E — Ylk=i anc ^ the right value 
of a is found by solving numerically the self-consistency 
equation 



aN 



E i 

E+(p)<0 



E i- 

E-(p)<0 



(21) 



This procedure can be repeated for all the lattice con- 
figurations shown in Fig. [2 obtaining again upper bounds 
for the energy of the ground state. In addition, we can 
compute the value of A at which the phase transition from 
the ferromagnetic ground state to the chiral state hap- 
pens. This is given by the configuration at which a^O 
becomes favorable, which is the configuration in which 
some E±(p) become negative. 



C. Numerical results 

In addition to the mean field studies, we have diago- 
nalized exactly the Hamiltonian in Eq. (JSJ) using all four 
geometries and an increasing number of spins. In all 
cases we observe a phase transition from a ferromagnetic 
regime (A = 0) to a chiral one (|A| — > oo). The location 
of the phase transition depends on the particular model 
and agrees with the estimates developed above. From 
the analytics, for the ladder of type A and for the ring 
geometry we find Atran ~ 1.1, while for the two other 
configurations we have Atran ~ 1.7. As shown in Fig. [2l 
these values are close to the actual location of the jumps 
in the chirality that are obtained numerically. 

It is interesting to note that, in the chiral phase, the 
total chirality strongly depends on the geometry of the 
lattice. Since our Hamiltonian favors the same orienta- 
tion of chirality on all plaquettes, the type-A ladder ex- 
periences a particular kind of frustration, where no two 
adjacent plaquettes are able to host maximal chiral cur- 
rents. A similar phenomenon is experienced by the type- 
B lattice, but not by the type-C ladder, where the weakly 
connected plaquettes are able to saturate the maximum 
value of the chirality per site (cf. Fig. [2]). On the other 
hand, the ring configuration appears to have only bipar- 
tite entanglement. This is due to rotational symmetry, 
which forces the central spin to decouple from each of the 
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FIG. 2: Absolute value of the chirality vs. strength of chiral 
interaction for the different geometries of Fig. [T] (for N = 9 
spins with periodic boundary conditions). We see that the 
chirality of lattice types A, B, C, and D becomes nonzero at 
the points Aa ~ 1.1, Ab ~ 1.7, Ac ~ 1.7, and Ad ~ 1.1, 
respectively. The horizontal lines (from top to bottom) cor- 
respond to the maximum chirality (= 2v3), and the entan- 
glement witnesses for W-like (= 3y/3/2), GHZ-like (= 2) and 
bipartite entangled (= 1) states, as demonstrated in Section 
II. 
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FIG. 3: Ground state energy per lattice site vs. strength of 
chiral interaction for a type-A ladder (Fig.[]J,). The solid line 
is the analytical estimate computed using the fermionization 
in Sect. IIH"Bl while the dashed and dashed-dotted lines come 
from an exact diagonalization of the model for N = 7 and 
N — 15 spins, respectively. 



triangles, thus forbidding any possible tripartite correla- 
tions. 

As we have seen, chirality acts as an entanglement wit- 
ness for two- and three-partite entanglement. Following 
our previous analysis, and from Fig. [51 we conclude that 
plaquettes in the A and B ladders have genuine three- 
partite entanglement of GHZ-type. Ladder C, with its 



6 




-10 -5 5 10 
A 

FIG. 4: Spin-spin correlations Co, a [Eq. (|22|l ] for a periodic 
type-A ladder with N = 24 spins. The dashed-dot line cor- 
responds to A = 0.2 and the solid line to A = 100. In the 
chiral regime the correlations decay exponentially (on either 
direction of the periodic ladder). 



weakly bound triangles and no frustration, achieves W- 
type entanglement on every three neighboring spins. 

We have compared the energies per site provided by 
the numerical simulation with the analytical mean-field 
results. As shown in Fig. [31 the derivative of the energy 
per site experiences a discontinuity around the transition 
point. This feature is shown both in the numerics and in 
the mean field estimates. The latter, however, only pro- 
vide a variational upper bound that is approached by the 
numerical exact diagonalization for increasing number of 
spins. 

We have also investigated numerically the spin-spin 
correlations of the ground state of type-A ladder for var- 
ious ladder sizes and chiral couplings, A, using periodic 
boundary conditions. A suitable way to measure the 
quantum correlations is to evaluate the connected cor- 
relator 

= (ffi ■ ffi) - • fa) (22) 

between any two spins. Scaling studies show that two 
different behaviors corresponding to the two phases of 
the system are obtained. Below the phase transition, 
A < Atran ~ on the ferromagnetic regime, spins 
are completely aligned and the connected correlation be- 
comes zero. Above the phase transition, in the chiral 
phase A > Atran ~ 1.1 the quantum spin-spin correla- 
tions appear to decrease exponentially fast with the sep- 
aration of spins, \i— j\. Both regimes are shown in Fig. 31 
where we plot the correlator Co, a between the first and 
any other spin on a type-A ladder with N — 24 sites. 
For A = 0.2 the correlator Co, a is identically zero, while 
for A = 100 it becomes negligible beyond the sixth site. 
This demonstrates the absence of structure in the ground 
state of type-A ladder when it is in the chiral regime. 
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TABLE II: Properties of the ground state in the chiral phase, 
A 2> 1, for various lattice sizes. We show the magnetization, 
S z , the total spin, S, and the degeneracy of the ground state. 



IV. TWO-DIMENSIONAL MODEL 

In this section we consider the two dimensional trian- 
gular lattice and impose periodic boundary conditions in 
both directions. The lattice is represented by an array, 
as illustrated in Fig. with dimensions denoted as N 
times M. We have performed exact numerical diagonal- 
ization of the Hamiltonian for various sizes of the lattice 
and characterized both the ground state and its lowest 
excitations. 

The results are summarized in Table [TT] and in Fig. [5] 
Once more, we obtain two phases separated by a phase 
transition point at a position that approaches Atran ~ 1.1 
for large lattices. When the chiral coupling A is small 
(A < Atran) then the Heisenberg interaction is dominant 
and the system is in a ferromagnetic state. When the 
coupling is large (A > Atran) then the chiral term is dom- 
inant causing the ground state to have non-zero chirality 
(see Fig. [5ja). Connecting with our studies of the chi- 
rality as an entanglement witness, the larger degree of 
frustration in the toric setup prevents the system from 
reaching a large value of the chirality, which lays around 
1.5 in the limit |A| — > oo. This implies that the state 
belongs to the set of states which are at least two-party 
entangled. 

The existence of the quantum phase transition and 
even its location agree qualitatively with the features of 
the Type A ladder, whose geometry most closely resem- 
bles that of these toric structures. However, as we show 
below, there are some differences. The most notable one 
is in the angular momentum of the ground state. We 
have found that in the chiral phase the system tends to 
adopt the state with the lowest total angular momentum 
which is compatible with the number of spins. Thus, as 
shown in Table [III if the number of spins is even, which 
includes the case of Type A ladders, the ground state 
is a state with S = and S z = and has no degener- 
acy. The corresponding momentum of the ground state is 
zero. However, if the total number of spins is odd, the to- 
tal angular momentum must be fractional, having S = | 
and S z = ±i. In this case the ground state becomes 
four-fold degenerate distinguished by the z component 
of the spin and by the momentum which is either zero or 
7T in both directions. 
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FIG. 5: (a) Ground state mean chirality as a function of the 
coupling constant A for different lattice sizes, (b) Energy gap 
between the ground state and the first excited state. 

It is interesting to note that the excited states almost 
always have a larger angular momentum and are sepa- 
rated by a large energy gap, 0(|A|), from the ground 
state. This is illustrated in Fig. [5jb). This energy gap, 
which survives in the thermodynamic limit, is responsi- 
ble for a finite correlation length. This causes both the 
spin-spin and chiral-chiral correlations to die off expo- 
nentially fast (see Fig. [4] for the spin-spin correlations in 
the case of a type-A ladder). 

Remarkably, on some special cases, such as the 3x5 
lattice, the gap between the ground and excited states 
is much smaller (see Fig. [SJb)). The excited states then 
have similar values of total angular momenta and their 
chiralities differ only slightly — by 0.2 or less — from 
that of the ground state. We conjecture that the exis- 
tence of these states is supported by the frustration of 
the chirality in Eq. ^ where we enforce all plaquettes to 
have a similar orientation of the spin current. 

The chiral-chiral correlations are presented in Fig.[6]for 
a 4 x 4 lattice, where we show the connected correlation 

(XX) = (XijkXi mn ) — (Xijk)(Xi mn ) (23) 

between a reference plaquette (Imn) and any other one, 
(ijk). In Fig.[5Ja) the reference plaquette is given by the 



b) 
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FIG. 6: A geometrical configuration of the 4x4 periodic lat- 
tice. The position (n, m) of the site on the lattice is given 
in terms of the integers n,m = 1,...,4. A certain bond is 
depicted as (n, m) — (k, I) that connects the corresponding 
two neighboring sites, (a) Chiral correlation between the pla- 
quette with sites (2,2), (2,3) and (3,3) and the remaining 
ones. Note the difference between the value of (XX) for the 
same plaquette (7.042) and neighboring ones, (b) Dimer cor- 
relations between the bond (2, 2) — (2, 3) and the remaining 
bonds. The simulations were performed deep in the chiral 
regime |A| = 100. 



sites (2, 2), (2, 3) and (3, 3) of the 4 x 4 lattice. The value 
of (XX) between the reference site and itself is 7.042, 
while the connected correlation with respect to neighbor- 
ing triangles decreases dramatically. This indicates the 
absence of any structure in the chirality of the ground 
state. 

The structure of the ground state is a singlet when- 
ever the total number of spins is even. This by itself 
does not necessarily imply a complicated hidden order, 
since we could achieve the same value of the total angular 
momentum by packing the spins into singlets. Following 
Ref. [2(| we have computed the dimer order parameter, 
which is defined as the connected correlation of an op- 
erator dij = (1 — (f; L ■ <Xj)/4 that detects singlets. The 
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FIG. 7: The triangular configuration. Three-spin interaction 
terms appear between sites 1, 2 and 3 since tunnelling be- 
tween two of the sites, such as 1 and 2, can happen directly 
or through the third site. The latter case results into an ex- 
change interaction between 1 and 2 that is influenced by the 
state of the spin at site 3. 



dimer-dimer correlations are then defined as 

(di 3 dki) - (dij)(d k i) 



(dki)(l - (dij)) 



(24) 



where (kl) is the reference bond and (ij) any other bond. 
In our case the value of the dimer order D decreases 
rapidly as a function of the bond separation (see Fig. [5^), 
suggesting that this is not a dimer solid. 

Summarizing, the ground state of the 4x4 periodic 
lattice appears to have no structure with respect to spin- 
spin correlations or chirality-chirality correlations. More- 
over, the dimer correlations appears to reduce dramati- 
cally for increasing distance between the dimers. 



V. OPTICAL LATTICE IMPLEMENTATION 

Here we would like to sketch a method of how to pro- 
duce the Heisenberg and chiral interactions studied in 
this paper by employing cold atoms superposed with op- 
tical lattices. More precisely we consider a Mott insu- 
lator of two species bosonic atoms loaded on an optical 
lattice with a triangular configuration, as shown in Fig. [7] 
These two states can be encoded in the internal hypcrfinc 
states of an atom. In the limit of deep off-resonance opti- 
cal lattices, the evolution of this system is described by a 
Bose-Hubbard Hamiltonian for atoms of species S = a.b 

H = -^(jfat s a i+ i, s + H.c.) 

+ - ^2 Us,srii,s(ni,s - 1) + U a ,bni, a ni t b(25) 



Here J s > is the atom tunnelling coupling between 
neighboring sites, Us,S' is the on-site interaction between 
atoms S and S", a^s denotes the annihilation operator 
of atoms of species S at the site i and n i: s is the corre- 
sponding number operator. For simplicity we take Us,s' 
to be positive. 

As we are interested in a chiral interaction we want 
to generate an effective charge-magnetic field coupling in 



the atomic system manifested by complex, position de- 
pendent, tunnelling couplings. For example, it has been 
shown in Refs. [lcil I 111. fl3 | that a gradient of a magnetic 
field along the direction of the electric dipole moment 
can simulate effectively the interaction between a charge 
and a homogenous magnetic field. Alternatively, Raman 
assisted tunnelling with lasers that have a phase differ- 
ence [13j can be also employed. We will employ such a 
technique here to generate, in a controlled way, complex 
position dependent tunnelling couplings. 

In the insulator regime, hopping is weak compared to 
the interaction, J -C U, and atoms have a low probabil- 
ity of jumping to other sites. For populations of only one 
atom per lattice site we can employ the pseudo-spin basis 
for each lattice site, given by \n a = 1, n& = 0) = |f) and 
\n a = 0, rib = 1) = II)- In that limit one can treat the 
hopping terms perturbatively and find that they help to 
establish effective exchange interactions between atoms 
in neighboring sites. Up to third order in the perturba- 
tion expansion the effective Hamiltonian is given by 

i i 
i 

+ effi ■ <7 i+ i x °i+2 (26) 
The presented couplings are given by 
J a2 J' 



J" 2 


(3 = 


ja2 


Ubb 


2U aa 


J" 2 


5 = 


(J a2 J l 


Ubb' 




1 \ m a 



'hb 



J al J b J b J a 



in 2 in 2 

Lu aa LtJ bb 



If we choose J b to be a negative imaginary number then 
it is possible to set a, 7 and S equal to zero. In addition, 
we can apply an effective magnetic field that cancels the 
/? term, resulting eventually to an isolated e term, 



(<?,+! X tTj+a). 



(27) 



This is the chiral interaction that we required. 

The presented interaction is produced from third order 
perturbation theory and it could be rather small [2l| . A 
similar interaction can be engineered using recently de- 
veloped techniques that involve cold molecules in optical 
lattices [l4| • Compared to the cold atom implementation 
shown here, the molecules would have the advantage of 
producing stronger interactions. 



VI. CONCLUSIONS 

In this article we have studied a chiral spin system 
which is defined on a two dimensional triangular lattice. 
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The analytical or numerical study of the 2D-model is 
hard so we initially resorted to simplified geometries such 
as triangular ladders and polygons (Fig.QJ. By fermion- 
izing and using mean field approximation we obtained the 
transition points, that separate between a spin-ordered 
and a chiral phase, as well as an upper bound for the 
ground state energies. These are in agreement with ex- 
act diagonalization of finite size systems (see Figs. [21 El- 
From these results we have observed that for the ladder 
of type A the phase transition happens for smaller values 
of the chiral coupling A, compared to the corresponding 
transition points for cases B and C. This has been at- 
tributed to frustration effects that also arise in the two 
dimensional case. On the contrary, case C has chirality 
that saturates the upper bound \ < 2\/3, for large chiral 
coupling A. Moreover, case A and the 2D-lattice saturate 
to a lower value of the chirality. Using the chirality as 
a witness of entanglement, we have discussed the entan- 
glement properties of the system and distinguished be- 
tween two-party and genuine three-party entanglement. 
We have suggested that the recent advances in cold atom 
and polar molecule technologies [l(| EH, G2, EH E3 could 
facilitate the experimental realization of the chiral sys- 
tems presented here. 

The study of topological order is rather complex and 
its presence in actual physical systems is hard to probe in 
a definitive way. Moreover, finite size simulations make it 



possible to deduce certain characteristics which can only 
indicate the presence of topological order. Indeed, the 
extensive numerical calculations carried out here on the 
2D model suggest the presence of topological order in the 
ground state of the system, which has been supported by 
a number of observations. Firstly, the system exhibits a 
ground state degeneracy that goes beyond the breaking 
of the PT-symmetry. Secondly, there is an energy gap be- 
tween the ground and the excited states, which persists 
with increasing lattice size. Finally, in the chiral regime 
the two-point spin correlations, the chiral connected cor- 
relation and the dimer correlations appear to decay ex- 
ponentially, a behavior that indicates the absence of any 
local structure in the ground states. 
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